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Abstract 

We use a characterization of the fractional Laplacian as a Dirichlet to Neumann operator 
for an appropriate differential equation to study its obstacle problem. We write an equivalent 
characterization as a thin obstacle problem. In this way we are able to apply local type 
arguments to obtain sharp regularity estimates for the solution and study the regularity of 
the free boundary. 

1 Introduction 

Constrained variational problems with fractional diffusion appear in the study of the quasi- 
geostrophic flow model ([9] and 0), anomalous diffusion 3J and American options with jump 
processes [10] (See also [16]). 

In this paper we will study the fractional obstacle problem. It appears in several contexts and 
it can be stated in different ways. Given a smooth function (p : R" M, that we can assume it 
decays rapidly at infinity, we look for a function u satisfying 

u{x) > (f for a; G R", 
{~Ayu{x) = for u{x) > ip, (1.1) 
{-Ayu{x) > for a; e R". 

This system of inequalities may also be written as a variational problem in (The Hilbert 
space spanned by functions with the norm ||w||^s = II^'(C)ICNIl2). The function u minimizes 
its norm in among all functions v satisfying v > ip. 

From a potential theoretic point of view, u can also be characterized as the smallest super- 
solution of {—Ayv > 0, among those functions v > (p. This makes sense since the minimum of 
supersolutions is also a supersolution. 

Finally, as a Hamilton-Jacobi equation we can describe u by the property 

min((— A)''u, u — ip) — 

Each of the previous descriptions is motivated by a different application and different inter- 
pretations suggest how to treat different issues as we study them. 

In this paper we develop two aspects of the problem: optimal regularity of the solution and 
regularity of the free boundary. The existing quasi-optimal regularity of the solution u was proven 
by one of the authors in [IB]. There it is shown that u E C^'" for every a £ (0, s) using methods 
mainly from potential analysis. 



In the particular case s = 1/2 it is easy to see that the operator (— A)^/^ coincides with 
the Dirichlet to Neumann operator in the upper half space of IR"+^. More precisely, given u{x) 
defined in R", extend it to u*{x, y) in M"+^ by convolving with the classical Poisson kernel. Then 
(-A)i/2y(x) =-u;(x,0). 

In [7], we obtain the same interpretation for any fractional power {—AYu (0 < s < 1) as the 
Dirichlet to Neumann operator of an appropriate extension u*{x,y). In this paper we use this 
characterization to obtain the sharp optimal regularity u S C^ '^. We also obtain a regularity 
result for the free boundary away from singular points. 

We obtain the equivalent problem to p.ip extending u to the upper half space u : M" x [0, oo) 

M. 

u{x,0)>ip{x) forxeM" (1.2) 

div (y°Vu(a;, y)) =0 for y > (1.3) 

hm y°'dyu{x,y) — for u{x,0) > Lp{x) (1.4) 

lim y^dyuix, y) < for x G K" (1.5) 

where s = (1 — a)/2. 

For y > 0, u{x, y) is smooth, thus the equation (|1.3p is understood in the classical sense. 
The equations at the boundary (|1.4|) and (|1.5p should be understood in the weak sense as it is 
explained in [7]. Also, from [16j . u{x, 0) G C^'°' for any a < s, in particular for some a such that 
2s<l + a<l + s, then limy^o+ y'^dyu(x, y) can be understood in the classical sense too. 

From |16| we know that u{x, 0) G C^'" for every a < s, which easily implies the same for 
u{x,yo) uniformly with respect to any fixed yo. We also know from '16] that for C = sup |Z)^(p|, 
deeU > — C for any unit vector e (m is semiconvex in the x variable). 

The function u that solves the problem above can be extended to the whole space by sym- 
metrization u{x,y) = u{x, —y). In 7 it was shown that the condition (jl.4p is equivalent to the 
symmetric extension of u being a solution of div(|y|°Vu(X)) — across the boundary {y = 0} 
on the part where u{x,y) > ifi^x). On the other hand, the condition (|1.5p is equivalent to 
div{\y\°-Wu{X)) < as a distribution. Let us call L^u = div(|y|°Vu). The setting for the 
symmetric extension translates as 

u{x,0)>ipix) fora;GM" (1.6) 

u(x,y) ^u{x,-y) (1.7) 

Lau{X)^0 ioT X eW^+^\{{x,0):u{x,0)^ip{x)} (1.8) 

Lau{X) < for X G M"+^ in the distributional sense (1.9) 

Notice that L^u is a singular measure supported on the set A := {{x, 0) : u{x, 0) = ip{x)}. By 
the continuity of u, we have that {u{X) — (p{x))Lau{X) = in R"+^. 

In the case s = 1/2 (i.e. a — 0), LaU — Au, and we are on the situation of the classical thin 
obstacle problem. Optimal regularity estimates, as well as free boundary regularity results have 
been obtained recently for this problem in [1], [2]. However, the results in this paper are new even 
to that case because we can consider nonzero obstacles tp. 

One of the great advantages of stating the problem with a PDF (|1.2m.5p instead of a nonlocal 
equation (jl.ip is that we can localize the problem. We will consider a local version of p.2m.5p . 
We write X = (x, y), with X G R"+\ a; G M", y G R, and 

Br = {X:Xe R"+\ |A| < r} 

B* = G R" : l^l < r} 

S"^ = {A : A G R"+MA| = r} 
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Given 95 : ^ K, we consider a function u : Bi ^ M. satisfying the following equation 



u{x, 0) > 1^9(2;) 

u{x,y) = u{x,-y) 
Lau{X) = 
Lau{X) < 



for X e Bi\{{x,0) : u{x, 0) = Lp{x)} 
for X e i?i in the distributional sense 



for X E Bl 



(1.10) 
(1.11) 
(1.12) 
(1.13) 



Compared to (II. 6111. 9|) . we are forgeting everthing that happens outside Bi. We can do that 
because we are dealing with a local PDE, and all the information that in the original problem 
(|l.ip would be coming from outside the unit ball is encoded in the values of u in Bi for y > 0. To 
study the regularity of the problem, we will focus on a solution to (|1. 10111. l3|) . 

Notice that the result of [1] is the particular case of this problem when s = 1/2 and ip = 0. 
The case (p ^ does not seem to follow from the case 1^ = in a straight forward way. The 
main idea of writing the problem as a local-type partial differential equation is to be able to use 
usual techniques for regularity of partial differential equations, like monotonicity formulas and 
classification of blowup profiles. 

The problem can be thought as the minimization of the functional 



from all functions v in the weighted Sobolev space W^'^{Bi, |j/|'^) such that v{x,0) > (p{0)- Fol- 
lowing the intuition explained in 7J this could be interpreted as an obstacle problem, where the 
obstacle is only defined in a set of codimension 1 + a, and a is not necessarily an integer number. 

The theory of degenerate elliptic equations in weighted Sobolev spaces plays an important role 
in this work, specially the results in the paper |12| . 

The strategy of our proof is the following. From , it is enough to show the regularity around 
free boundary points, so we assume that the origin (0, 0) belongs to the free boundary. We will 
show that v?' lyl^dtr decays in the appropriate way. To see that, we will study the limit as 
r — > of a variation of the Almgren's formula: 



In order to obtain the possible limits of Almgren's formula, we make a blowup and it turns 
out that the limits correspond to the degrees of the global homogeneous solutions of (|1.6m.9[) with 
= 0. Finally. The minimum possible of such degrees of homogeneity can be found either using 
a monotonicity formula similar to the one from [1] , or using the sharp result of [16j for the case of 
convex contact sets. 

Later, the regularity of the free boundary is addressed using that in the blowup solutions the 
contact set is a half space. 

2 Preliminaries 

We will start with some elementary properties of the equation LaV = 0. In the case a = 
most of these properties are very classical results for harmonic functions. Intuitively, these results 
generalize to the case a G (—1, 1) through the idea of the fractional-dimension extension [7]. On 
a first reading, it may be a good idea to skim through this section quickly. 

Since |y|° is an A2 weight, the following proposition is a particular case of a result in T2] that 
we will use in this paper 




$(r) 
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Proposition 2.1. Assume LaV = in Br- Then there is an a > such that the function v is 
C°'{Br/2) and 

„ n C 

For a proof of the above proposition check [12] . 
The Harnack inequality is also available. 

Proposition 2.2. Assume LaV = in Br, then 

sup u < C inf u 

Br/2 

Moreover, if instead we had a right hand side: LaV — f{X) in Br, we would have 

sup u < C inf u + Cr^ sup I/I (2.1) 

The estimate with the right hand side (|2.ip can be deduced from the weak-Harnack inequalities 
apphed to u + c|Xp and u + c{r'^ — |Xp), as it is standard. Notice that this construction can be 
translated in the x direction but not in y. We are goint to apply the estimate (|2.ip only in balls 
centered on {y = 0} (in the proof of Corollarv l7.4p . 

Using Proposition 12.11 and the translation invariance of the equation in the x variable, we 
obtain the following result. 

Proposition 2.3. Assume LaV = m BriXo) for some r > 0. Then 

sup \VxV\ < — OSC V 

Br/2(Xo) r BAXo) 

c 

[^xv]c''(B,,AXo)) < ^osc V for some small a>0 

where write [/]c"(-D) to denote the seminorm 

\f{^)-f{y)\ 



[/]c°(D) = sup 



,yeD \x-y\°' 



Remark 2.4. Notice that the estimate in Proposition 12.31 refers to the derivatives with respect 
to X, not with respect to y. 



Proof. Assume first r = 1. From Proposition 12.11 we know apriori that v is C" and 

II^IUafR fy \\ ^ C OSC V 

Given a tangential unit vector r = (t',0) The incremental quotient: 

v{X + Ht) - v{X) 



Vh 



\h\^ 



is bounded in B-j/s independently of h. Moreover, by the linearity and translation invariance of 
the equation, Vh is also a solution of LaVh = and then by Proposition 12.11 it is C". 

This implies that v is of class C^" in the x variable (see [5] , or [T3] to see an example where the 
method is only applied in the tangential direction). Iterating the same argument a finite number 
of times we obtain that v is Lipschitz and finally C^'" in the tangential directions in -B1/2 

The estimate of the Proposition follows by scaling. □ 
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Corollary 2.5. Assume LaV = in Br^Xo) for some r > 0. Then 

C 

sup \D'1v\ < —r osc V 
C 

\Dlv]c''{B^/,iXo)) < -j^ osc V for some small a>0 
Proof. Assume first r — I, the general case follows by scaling. By Proposition 



sup \'^xv\ < C osc V 

Moreover, Wj — dx^v is also a solution to LaWj = in Bi{Xo) since the equation is linear and 
translation invariant in x. Applying Proposition 12 . 31 again we obtain 

sup \DIv\^ sup ^\^xWj\ 
Bi/i(Xa) Si/4(Xo) J 

< C osc w 

Bi/2{Xo) 

< C osc V 

Bi{Xo) 

For k — 2, the proposition follows by scaling. 

Iterating Proposition 12.31 further, we obtain the bounds for larger values of k with the same 
reasoning. □ 

Proposition 2.6. Assume LaV = in Br{Xo). Then for any r < I, 



sup 

Br/2{Xo) 



a 

y 



c 

< osc V 
H baxq) 



Proof. From Corollary 12. 51 we have that 

C 

sup \Dlv\ < — osc 



However, from the equation we have that 



A 

L\x'^ — "^yy "^y 



Therefore 

I a . C 

sup \Vyy + -Vy \ < ^ OSC V 

B^M^o) y r^ BAXo) 



The following is a Liouville type result 
Lemma 2.7. Let v be a global solution of 

Lav{X) = X e K" X M 
such that v{x,y) — v{x, —y) and \v{X)\ < CjA"!*^. Then v is a polynomial. 



□ 
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Proof. We use Proposition 12.31 and Proposition 12.61 and induction in the degree k. The following 
elementary fact will be used: if Vj;U is a polynomial and v{Q,y) is a polynomial in y, then w is a 
polynomial. 

Let us start with the case A: < 1. By taking r — > +oo in Proposition l2.6l we obtain Vyy + ^Vy = 0. 
This gives a simple second order ODE for each x whose solutions have the general form by\y\~"' + c. 
Therefore, for a fixed x, v(x,y) — by\y\~'^ + c, but since we assume that v is symmetric in y, we 
have that v{x,y) must be a constant for every fixed x. 

On the other hand, taking r +oo in CoroUarv 12.51 (with fc = 2) we have that D'^v — 0. 
Therefore for each fixed y, v{x,y) is a first order polynomial in x. 

Combining the two facts above we have that if A: < 1 then d is a polynomial of the form 
v{x, y) = bx + c. 

Now we consider larger values of k. 

From Proposition 12. 3[ |Va;i'(X)| < ClXj*^^^. Moreover, is also a global solution of the 
same equation symmetric in y. By the inductive hypothesis is a polynomial of degree k — 1. 
From Proposition [2Jl \vyy + ^Vy\ < C\X\''^^. Observe that 

Vyy + \y = | J/ " (M^dyV) 

Thus Vyy + ^Vy satlsflcs the same equation as v {\y\'^Vy satisfies the conjugate equation. See [7])- 
Therefore Vyy + ^Vy is a polynomial of degree k — 2. 

In particular, for a; = 0, Vyy{0, y) + ^Vy{0, y) is some even polynomial p{y) of degree at most 
fc — 2. Let us say that p{y) = oq + 02?/^ + • • • + a2dy'^'^- Then u(0, y) must be 

v(0 v) = c + bvlvr" + - + - + ■■■ + — 

v[U,y) c + Dy\y\ + 2{l + af ^ A{3 + af ^ ^ (2d + 2)(2d + 1 + a) ^ 

but since v is even in y, b = and v{0, y) is a polynomial. 

Since ^(O,?;) is a polynomial in y and V xV is a polynomial, then w is a polynomial. □ 

Remark 2.8. The symmetry condition v{x,y) — v{x,—y) is necessary. The simplest counterex- 
ample without that condition is v{x,y) = \y\~'^y. 

The equation LaU — can be understood as the Laplace equation in n + 1 + a dimensions [7]. 
In that context, the following lemma is just the mean value theorem. 

Lemma 2.9. Let v be a function for which 

Lav{X) < C\yr\Xf forXeBi 

Then 

«(0) > / viX) ly^da - Cr^+^ 

for any r < 1. where 



i-'n+a 



Si 

Proof. The proof is essentially the same as the mean value property for harmonic functions. Let 
us consider the case C = first. 
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We consider the test function 

T^f v\ ( Cn,a Cn a „\ 

r(A) = max 7——— — r- — -,0 

where Cn,a — + a — l)^'^uj^^]^g^. Note that F is supported in Br, F > and LaT{X) ~ —6q + /i, 
where jjL is the measure supported on Sr given by r~"~^w^]^g\y\°' dcr (See [7]). Thus 

< j -Lav{X)V{X) dX 
= jvv{X)-VT{X) \y\^dX 

Bi 

= J -v{X)LaT{X) dX 

Bi 

= «(0) I «(^) l^^l^d'^ 

Sr 

which proves the Lemma when C — Q. In the case C > 0, we apply the above computation to 
v{X)-C\X\'^/{2{n + a + l)). □ 

We will also need Poincare type inequalities in the context of weighted Sobolev spaces. The 
following is a classical Poincare inequality whose proof can be found in [12j . 

Lemma 2.10. For any function v G W^'^{Bi, \y\°') the following inequality holds 

j |«(X)-iJ|2 |yrda<Cr j \Vv{X)\' \y\'^dX 

Sr Br 

where v ~ ^ ^ \.n+a /g "^^i^) Ij/I^dcr and C is a constant that depends only on a and dimension. 

Remark 2.11. In [T^], the inequality is done with v being the average in the whole ball instead 
of the sphere. Indeed, the modification in the proof is straight forward given that a function v in 
the weighted Sobolev space W^'^{Bi, jyl") has a trace in L'^{Si, \y\'^)- 



Another form of Poincare inequality that will come handy is the following. 

Lemma 2.12. For any r < 1 there is a constant C > (depending only on r, a and dimension) 
such that given any function v € W^'^(Bi, \y\°') the following inequality holds 

j \v{X) - v{rX)\'' |yrda<Cr j \Vv{Xf \y\''dX 

Si Bi 

The proof is standard. It can be done for example integrating Vw along the lines sX with 
s € (r, 1), or using the compactness of the trace operators from W^'^{Bi, to L'^{Sr, 

Since the main difficulties in our problem appear in the behaviour of the solution near the free 
boundary, we will assume that the origin belongs to it. 

Almgren's monotonicity formula was proved for the problem 10111.1^ in [7\ if (p — 0. We 
cannot reduce the problem to that case. Instead, assuming (p G C^'^, we let u(x,y) = u{x,y) — 
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ip{x) + -^j^^y^ so that LaU = at the origin. Denote hy A = {u = 0} = {u = ip} the contact set. 
The function u has the following properties: 

u{x, 0) > for X £ Bl (2.2) 

u{x,y) = u{x,-y) (2.3) 

Lau{=c, y) = |yr (A^(a:) - ^^{0)) =: \y\^g{x) for {x, y)eBi\A (2.4) 

LaH^, y) < ly^gix) for {x, y) e Bi (2.5) 

where we define the function g{x) = {A(p{x) — A(p{0)) that is Lipschitz as long as </? € C^'^. 

The only problem we have is that the right hand side in (12. 4|) is not zero. However, \y\"'g{x) = 
\y\°'{Aip{x) — A(p{0)) that decays to as a; ^ 0. In fact 

\Lau{x,y)\<C\ynx\ (2.6) 

We would expect a small variation of Almgren's monotonicity formula to work, or at least to 
remain bounded as r — > when we apply it to w. The functions u and u have the same regularity 
respect to x, so we will prove our estimates in terms of u. In order to simplify the notation we 
will write u from now on, but we mean u. 

Lemma 2.13. Let u be a function for which ((^ . ([^ . (^3)) hold and u{0) = 0. Then 

J \u{X)\^ ly^daKCr J \Vu{X)\^ lyl^dX + Cr^+''+" 

Sr Br 

for a constant C depending only on a, n and \\(p\\ij2,i- 

Remark 2.14. The exponent in is a consequence of the chosen regularity for tf, in this 

case C^'^. 

Proof. By Lemma |2.10[ we have 

J \uiX)f \y\-da <CrJ \Vu{X)\' ly^dX + u J u{X) \y\''da (2.7) 

Sr Br Sr 

Therefore, to prove the lemma we need to find a suitable upper bound for 

^ \u{X)\ \y\^da 

Sr 

From Lemma 12.91 we have that 



= u(0) > ^ / u{X) \yrda - Cr^ 



therefore 

' u+{X) |j/|°d(7 < / u-{X) |y|Mcr + Cr3+'^+" (2.8) 



Sr 
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Now we estimate u {x,y) by integrating along the straight hne (2^,0) to {x,y) and applying 
Cauchy-Schwarz. Since u{x, 0) > 0, 

y 



/ y \ / y \ 1/2 

<{ j \Vu{x,t)\^t'' dt\ iji-^dt 

y \ 1/2 



< Cy^ J |Vu(x,i)|^i'' dt 



Now we integrate the above inequality on 5*^ and apply Cauchy-Schwarz again. Notice that on 
5*^ we have da = "^dx. 



u'{x,y)\y\°'da ^ j u {x,y)\y\'' r dx 

y 



1/2 



<Cr I |y|"T- I j |VM(x,t)|^r dtj dx 

1/2 / ^ 1/2 

<Cr\ I \Vuix,y)\' ly^dX | ( / \yrUx 

1/2 

< C'^(n+i+a)/2 [ I \yu{x,y)f jy^dX 
Combining with p.Sp we obtain 

1/2 

\u{X)\ |2/|MfT < Cr("+i+'^)/2 I / \yu{x,y)f jy^dX I + Cr3+'^+" 



Sr 

Putting this estimate back in (|2.7p we have 

\u{X)\-' \yrda<Cr I \Wu{Xr Ivl^dX + ^^ | / \u{X)\ {y^da 

<Cr j \Vu{X)\^ |y|"dX + Cr6+''+" 

Br 

□ 

Corollary 2.15. Let u he a function for which ([^ . ([^ . ((23), ^old and w(0) = 0. Then 

\u{X)\^ \y\"-dX<Cr'' j \Vu{X)\' jy^dX + Cr^+"+" 

Br Br 



for a constant C depending only on a, n and \\ip\\^2,i-- 

Proof. Take Lemma 12.131 and integrate in r. □ 

3 Frequency formula 

A crucial ingredient for our blowup analysis is the monotonicity of a frequency formula of Almgren 
type for functions u that solve (|2.2II2.5[) . and u(0) = 0. In the special case that the right hand side 
in ()2.4|) is zero, then the frequency formula takes the simple classical form as shown in [7j: 

r/^ \Vuf \y\-dX 
r) — ^ 

J,^ \u\' \y\^da 

In order to account for the right hand side in ()2.4p . we have to modify the formula with a 
suitable lower order term. We define the function: 

Fu{r) = J u{X)' \yrda (3.1) 

Notice that in terms of F we have 

r— logOF„ r = , , ,2 , , , +n + a. 

Thus, the classical frequency formula, shown in 7 , can be rephrased as that r-^ log(i^„(r)) is 
monotone nondecreasing in r. We will use the following modification 

= (r + Cor2)-ilogmax(F„(r),r"+'^+4) (3.2) 
dr 

As long as we are talking about only one function, we will write $ = Whenever there 
might be some ambiguity we will use the subindices. 

Theorem 3.1. For some small and a constant Cq large enough (depending on a, n and 
\\f\\Q2,i), The function $ (given in \S.2\) is monotone nondecreasing for r < r^. 

This result is an extension of the corresponding one in [7 . The proof is somewhat technical 
and we believe it is better to skip it on a first reading of this paper. In order not to distract the 
reader we postpone the proof to the appendix, so that we can quickly concentrate on the regularity 
of the obstacle problem in the following section. 

4 Local C^'" estimates 

In 'Iff, it is proved that the solutions of (II. ip are in general C^'" for any a < s. We can take this 
result to the context of ()1. 10111. 13|) locallizing the problem using a cutoff function. 

Lemma 4.1. Given a function ip : B'l — > M, we consider a function u : i?i — > M satisfying 
equations li2.S^2.5\) . Assume 

\uiX)\<C forXeBi 
\ip{x)\ + 115^(^)1 <C forxe B* 

Then we have the following estimates depending only on C , a and n: 
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• Uf,e{x, 0) is hounded below for x G B*^/2 f'^'"' '^^V vector e such that e ■ e„+i = 0. 

• lim;,^o \y\''Uy{x,y) G C"{Bl^^) for any a <1- s. 

• u{x, 0) £ C^'°'{B^^2) for any a < s. 

Proof. We apply here the results in [IB] . The only technical difficulty in order to apply immediately 
the main theorem therein is that our function u has a bounded domain instead of being defined in 
the whole space. The solution is to use a smooth radially symmetric cutoff function 77 : R"+^ M. 
that vanishes outside Bi and equals 1 in -62/3- 

We want to show that u{x) = ri{x,0)u{x,0) satisfies the hypothesis of Proposition 5.7 in [16]. 

Since the cutoff function 77 is radially symmetric, rju keeps the same Neumann condition on 
the boundary R" x {0}: 

lim y"'dy{riu){x,y) — for u{x,Q) > 

lim y''dy{r]u)ix,y) < for a; e K" 

However, the function rju is not necessarily a solution of Lafju = in the half space R" x (0, +00). 
Let f{X) — La {ri{X)u{X)). We solve the following "correction" equation: 

w{x,0)^0 forxeM" 
Lawix, y) = f{x, y) for y > 

We look at the right hand side f{X) — La {r]{X)u{X)) = wL^Ty + y'^VuVr]. Notice that 
f{X) = for X G -Bi/2. From Corollary 12.51 we have that w is (B*^^), and thus also is 
liiTiy^o \y\°'dyw{x,y). 

Now the function rju — w is a solution to 

Laivu - w){X) = 
in the whole half space. Its restriction u{x) = ri{x, 0)u{x, 0) — w{x, 0) satisfies 
u>0 

{~Ayu{x) > 4>{x) := -c lim y''dyw{x,y) 

{—A)'^u{x) = 0(a;) for those x where > 

Moreover, for (p = (— A)^'*^, u + ip is a solution of p.ip . In P2] it is shown that (— A)''u 
is bounded and Uee is bounded below for any unit vector e. From the first one we obtain the 
boundedncss of liuiy^o \y\°'Uy{x, y) and from the second one the bound below for Uee{x^ 0). 

Applying Proposition 5.7 in [16] . we get u £ C^ " for any a < s in the x variable and ob- 
tain the other estimate. Recall that (— A)'*u(a;) = limj,^o+ y°'9yu{x,y), so the C" estimate for 
lim^^o \y\°'Uy{x,y) comes from the corresponding estimate for (— A)*u(a;) in |16j . □ 

Remark 4.2. Notice that the function w that we construct in the proof above complies with all 
the hypothesis (4.1-4.3) and (4.32 - 4.34) from [16] that are needed in order to apply Proposition 
5.7 from there. In |16| only obstacles ip with compact support are considered for simplicity in the 
proof of existence of a solution, but the results extend whenever the problem makes sense. 

A C^'" regularity estimate with respect to the x variable only implies the same regularity 
estimates in both x and y for the case a = 0. In any case, what we can say is that y°'dyu{x,y) 
remains uniformly bounded as y ^ 0. 
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Proposition 4.3. Let u be as in Lemma \4.1\ Then we have the following estimates 



V,u(X) e C"(Bi/2) for any a <s (4.1) 

Ivl^dyuiX) € C^iB^f^) foranya<l-s (4.2) 

and the estimates on the corresponding norms depend only on the constants in the statement of 
Lemma \4-.1\ 

Proof. From Lemma l4.1l we have both estimates at the boundary {y = 0}. Since Wi — dx^u also 
solves LaW = in Si \ {?/ = 0}, the estimate extends to the interior and we obtain (|4.ip . On 
the other hand w{x.y) — \y\°- dyu{x ^ y) solves the equation div(|j/|~''Vw(X)) = (See jTj) and we 
obtain (g^l). □ 

The above proposition will be used together with the following compactness lemma to prove 
the existence of a blowup limit. 

Lemma 4.4. Let vj : Bi ^ M. be a bounded sequence of functions in Assume 
there is a constant C so that 

\yxVj{X)\<C forXeBi 
\dyV,{X)\<C\y\-'' forXeBi 

and also that for each S > the sequence is uniformly C^'" in Bi^g H {|y| > S}. 
Then there is a subsequence Vj^ that converges strongly in W'^''^{Bi, |y|°). 

Proof. Since the sequence Vj is bounded in W^'^{Bi, |y|°), then there is a subsequence that con- 
verges strongly in L^{Bi, \y\'^) (see [13]). 

Since for each 6 > the sequence is uniformly C^'" in Bi^s n {\y\ > S}, we can extract a 
subsequence so that Wvj converges uniformly in Bi^s ^ {\y\ ^ ^} for any S > 0. Thus, Wvj 
converges pointwise in Bi \ {y = 0} . 

In order to get strong convergence in W^''^{Bi,\y\°-), we wiU show that each partial dkVj 
converges in L^{Bi, \y\'^). Notice that each partial dkVj already converges almost everywhere 

We use two facts. The first is that if a sequence of functions Wj is bounded in L^ (for 1 < p < oo) 
and converges almost everywhere, then it converges weakly in L^. The second is that if Wj is 
an almost everywhere and weakly convergent sequence of functions bounded in LP{Q) for some 
measure space (ri,/i) such that /^(f^) < +oo and p > 2 then Wj converges strongly in L'^(n) (See 

Let us do it first for dyVj. Since \dyVj\ is bounded and converges in Bi\ {y = 0}, then 
\y\°'dyVj converges strongly in L^{Bi), and therefore Vj converges strongly in L^{Bi, jyl"). 

With respect to WxVj(X), we know that it is bounded uniformly in L°°{Bi, jj/l") and converges 
pointwise in _Bi \ {y = 0}, therefore it converges strongly in L^{Bi, |y|°). □ 

Remark 4.5. Notice that the function u that solves (|1.10m.l3|) can only be C^'" in both variables 
x,y in the case a < 0. If a > 0, then from the fact that y°'dyu{x, y) has a nonzero limit as y — > 0^ 
for some points x in the contact set {u = ip\, then dyU can never be bounded. However, after the 
change of variables z = y\y\~°' (see [7]), it is to be expected that the function u would be C^'" in 
X and z for a certain value of a. 
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5 Blowup profiles 



The purpose of this section is to characterize all posible blowup profiles of solutions to the equation. 
The problem we consider is to find all functions u : K"^^ B. with the properties: 

u is homogeneous of degree k (5-1) 

u(a;,0)>0 forxGR" (5.2) 

u(x,y) ^ u{x,-y) (5.3) 

LauiX) = in M" X M\ {(a;,0) : m(x,0) = 0} (5.4) 

Lau{X) < in M" X M in the distributional sense. (5.5) 

UrriX) > for any vector r in R" x {0} (5.6) 

The following proposition gives a lower bound for the degree of homogeneity fc, and it is all 
that is needed for the optimal regularity result. 

Proposition 5.1. // there is a solution u of lid. 1^5.6]) . then k > (3 — a)/2 = 1 + s, where 
(1 - a)/2 = s 

Proof. The idea of this proof is to use the main theorem from 16J. The theorem cannot be 
applied immediately because our function u is not bounded, and even worse, in general its growth 
at infinity is too high to even define the fractional Laplacian as a distribution. The solution is to 
use a smooth radially symmetric cutoff function r/ : M"+-'^ M. like 



1 if \X\ < 1/2 
-1^1' if IXI > 1 



Now we want to show that u{x) = r]{x,Q)u{x,0) satisfy the hypothesis of theorem 5.2 in 
[TB]. Notice that since u is convex in the x direction, then the contact set {(a;,0) : u{x,0) = 
0} = {(a;,0) : u{x,0) = 0} is convex. The proof follows as in Lemma [4. II but applying Theorem 
5.2 in [T^. Applying that theorem we conclude that u £ C^'^ in the x variable, and therefore 
fc > 1 + s. □ 

Remark 5.2. An alternative approach to show Proposition 15. II could be to use a monotonicity 
formula similar to how it is done in ^1^ . 



Even though the proposition above is enough to obtain the optimal regularity of the solutions u 
of (|1. 10111.13]) , we need to fully characterize all solutions of (|5.1II5.6|) in order to study the regularity 
of the free boundary. 

As before, we have A = {(a;, 0) S R" : u{x, 0) = 0}. Let us also define 

A, = {(x,0) £ M" : lim y''dyu{x,y) < 0} 

Notice that A* is the support of LaU, and by (|5.4p A* C A. Moreover, by the continuity of u, A 
is closed and A* c A. 

Lemma 5.3. // A* has H^^ -measure zero then u is a polynomial of degree k. 
Proof. From Proposition 14.31 \v\°'Uy{x. v) is locally bounded. If A* has _ff"-measure zero, then 

lim \y\°'Uy{x,y) — a.e. in x. 
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Therefore for limj,_+o \y\'^Uy{x, y) = weak-* in L°°. In [7], it is shown that this imphes that u is 
a global solution of 

Lau = in R" X R 

We conlude the proof using Lemma [27fl □ 

Proposition 5.4. If A has positive H^-measure then either u = Q, or k = 1 + s and A is half of 

W\ 

Proof. We first observe that if A* has i?"-measure zero then m = 0. Otherwise from Lemma 
15.31 uix. 0) would be a polynomial vanishing in a set of positive measure in R", thus it would be 
constant zero. The polynomial u must have the form 

u{x,y)^pi (x)y^ + ... Pi {x)y'^^ 

Computing LqU in terms of the above expresion, a simple iterative computation shows that 
LaU = only if pi = p2 = ■ ■ ■ = Pj = 0. 

Let us now consider the case when A^, has a positive i/"-measure. 

Since u is homogeneous, A* is a cone. Let us assume that e„ is a direction inside this cone. 
Since A* has a positive _ff"-measure, then it must be a thick convex cone in the sense that a 
neighborhood of e„ is contained in A*. Therefore, for any x G M", x + hen G A* for h large 
enough. 

Since 

u{x,y)-u{x,0) .. . .a , X ^ „ 
c lim p— = Imi \y\ Uy{x, y)<0 

y^a \y\ "-y y^O 

then u is not positive in a neighborhood of e„. Therefore, for any X E M"+^, u{X + he„) < for 
h large enough. 

By hypothesis (|5.6p . u is convex in the e„ direction. The function w ~ is decreasing and 

cannot be negative at any point X because otherwise lim/j^+oo u{X + hen) — +oo contradicting 
the above. 

On the other hand, w = on A and 

Law{X)=Q in R"+i \ a; D M"+i \ A 
Thus, w must be the first eighenfunction corresponding to minimizing the spherical integral 

\Vew\^\yrd0 



from all functions w such that i« = on A and 



Sr 



Since A is convex, it covers at most half of the sphere Sr Ci {y = 0} . If it was exactly half of 
the sphere then we have the explicit expresion 



■w{x, y) ^ c i^-s/xl + y2 - 



which is, up to a multiplicative factor, the only positive solution of LaW = which vanishes in 
{y = A x„ > 0} (the computation is somewhat lenghty, but it can be done quickly with a 
computer algebra system). 
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Notice that the above expUcit function is not a solution accross {y — /\Xn > 0}. Therefore for 
any convex cone A that is strictly contained in = A x„ > 0}, there must be another function 
which gives a smaller eigenvalue, and then a smaller degree of homogeneity k than 1 + s. But from 
Proposition 15.11 k > 1 + s. Therefore the only possibility is = 1 + s. Moreover, we also see in 
the explicit function w that A has to be {y = A x„ > 0}, which is half of R". Recall that e„ is 
an arbitrary direction inside the cone A. □ 

Proposition 5.5. Up to rotations and multiplicative constants, there is a unique solution of i5.1\ - 
\5.6\) that is homogeneous of degree 1 + s. 

For this solution the free boundary is fiat, there is a unit vector e such that A = {(a;, 0) : x • e > 

0}, and deU — c (x ■ e)^ + y^ — {x ■ e)j . 

Proof. In the proof of Proposition 15.41 once u^^ is uniquely determined as 



W = -Ux„ = C \^\/ xl + J/2 _ Xnj 

we have A = {u = 0} = {y = A x„ > 0}, and integrating on the lines parallel to e„ we determine 
u{x, 0) for every x. Now, if we had two solutions ui and U2 homogeneous of degree 1 + s, coinciding 
on {y = 0}, then necessarily ui(x,y) — U2{x,y) = c|?/|*y for some constant c and y > 0. But that 
constant must be zero in order for ui and U2 to be solutions accross {y = 0} \ A. 

Replacing e„ by an arbitrary unit vector e normal to dA we get the rest of the result. □ 

Corollary 5.6. Let u he the solution of J5. jH5.6)) that is homogeneous of degree I + s and such 
that e„ is normal to its free boundary, then u is constant with respect to Xi, X2, ■ ■ ■ ,Xn-i- 

Proof. Translating u in any direction r orthogonal to x„ and y would give another global profile 
with the same free boundary. By uniqueness, u must be invariant in those directions. □ 

We can summarize the results of the section in the following Theorem 

Theorem 5.7. Let u be a solution of i5.1Vl5.6]) . There are two posibilities: 

1. The degree fc = 1 + s. A is a half space, and u depends only on two directions. 

2. The degree k is an integer greater or equal to 2, u is a polynomial, and A has H" -measure 
zero. 



6 Blowup sequences and optimal regularity 

The optimal regularity of the solution will be obtained by carefully analyzing the possible values 
of $(0). Recall that when <i> is constant, its value depends on the degree of homogeneity of the 
function. In this case <&(0) represents the assymptotic behaviour at the origin. 

Lemma 6.1. Let u be a solution to li2.<M2.5\) . Then either ^^(O) = n + a + 2(1 + s) or $„(0) > 

n + a + 4 

In order to prove Lemma |6. II we will consider the following scaled version of F that represents 
the growth of u at the origin. 
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Using the value of dr we define the following nonhomogenous blowup sequence that will also 
be useful for studying the regularity of the free boundary in the next sections. 



Ur{X) = -^u{rX) 

dr 



(6.1) 



Lemma 6.2. Let u be as in Lemma \6.1[ //liminf,.^o % = +00, then there is a sequence —> 
and a nonzero function uq : M. such that 



uniformly in B1/2 
uniformly in -B1/2 
uniformly in -B1/2 



(6.2) 
(6.3) 
(6.4) 
(6.5) 



Moreover, uq satisfies the conditions id. 1^5.6]) and its degree of homogeneity is (<i>„(0)— n— a)/2. 

Proof. The function Ur is constructed in a way so that ||Mr|li2(5^ |j,|a) = 1 for every r. Since 
liminfr^o 75- = +00, then F„(r) > 7'"+"+4 fpj- j. small enough. Let us also consider the values of 
r smaller than the tq of Theorem 13.11 

We will first show that Ur remains bounded in W^'^{Bi, |?/|°) using the monotonicity of the 
frequency formula. We have that 

$(ro) > $(r) 

> (r + Cor2)^logmax(F„(r),r"+''+4) 
dr 

>{r + Cor^) + {n + a){l + Cor) (since F„(r) > r"+'^+4) (6-6) 



> 2(r + Cor^) 



+ {n + a){l + Cor) 



Recalling that \Lau\ < C|y|°|a:| in R"+i \ (A x {0}), and that w = on A x {0}, we estimate 

1/2 



uLaU dX 



< Cr""^ 



By Corollary [HTg 



\u{X)\^ \y\''dX<Cr^ / |Vu(A:)|^ jy^dX + Cr 



7+a4-n 



Thus 



uLaU dX 



1/2 



Ti + a + 5 



|Vm|2 \y\-dX 



Cr 



(6.7) 



From Lemma 12.131 



j Y^uiX)]^ \yrdX>^- J \uiX)\' ly^da^Cr 



,5+0'+'"' 
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but since we have dr/r^ — > +00, 



Sr 

for a constant C as large as we wish as r \ 0. Thus 

|Vii(X)|^ Ij/l-'dX > C'r"+''+3 - Cr'^+'^+" (6.8) 

Br 

> ^r"+''+3 for r small (6.9) 

Putting it all back in (|6?7|) . 



uLaU AX 



< Cr"+"+'* for r small 



Comparing with (|6.9p we see that for r small 



Br 



Now we continue with (|6.6p . 

cl>(.o) > (. + C^r^y-fj^T^ + + + ^or) 

Bi 

Combining the above inequality with the fact that /^^ |urP lyl^du — 1 and Poincare inequality, 
we obtain that the sequence Ur remains bounded in W^'^{Bi, \y\'^). 

Since = max.{ur,0) and u~ = max(— Wr,0) are subsolutions of the equation 

div(|yrV«) > -Cr\y\''\x\ 

then we have that Ur is bounded in L°°{B^/4) (See [H])- 

We also know from Lemma |4 . 1 1 that u is semiconvex in x, or in other words that d^eU > —C 
for any tangential unit vector e. This implies a bound of the same type for all the sequence 
since we have the choice of scaling lim^—^o 75- = 00. 

Notice that the functions Ur are solutions of a uniformly elliptic equation with smooth coeffi- 
cients in BiD {y > 6/2} for any 5 > 0. Then clearly u,. are uniformly C-^'" in Bi^s H {y > (5} for 
any a,6 > 0. We can therefore apply Proposition 14. 31 and Lemma [44l to obtain a subsequence u,.. 
such that Ur^ converges strongly in W^^'^(i?i/2, lyl") to some function uq. 

From Theorem 13. H we have that $u(f") is monotone and converges to some value as r — > 0. 
Moreover, the computations right above show that as r — > 

$„(rs w rs ^ + {n + a) 

~ — f — 1 — nil, \- [n + a) 
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Now we let s = rj ^ 0, to obtain 



$„(0)-(n + a) (6.11) 



In order to pass to the limit in the above expression we need Ug — > uq in W^'^{Br, |y|°) (which 
we have), and we also need that the denominator remains bounded away from zero (i.e. that 
uo 7^ 0). When r > 1 — (5, if (5 is small, this is a consequence of Lemma [2.121 In other words, we 
have that if 6 is small enough 



\Us\' 



for a constant c depending on 6, $(1) and dimension. 

We can iterate this inequality k times so that {1 — 6)'^ < r and obtain a uniform bound from 
below for the denominator (|6.10p . so that we can pass to a nonzero limit and get (|6.1ip . 

Since limr_>o ^ = +oo, we have that for any unit vector r = (t', 0) 

drrUr = —UrrirX) > -~C— > 

So, in the limit, uq is convex in the x direction. 
Each function Ur is a solution to 

Ur{x,0)>0 for a; e Si (6.12) 

LaUr{X) = '^Lau{rX) = ^\y\''girx) for X G Si \ {(x, 0) : u,(a;, 0) = 0} (6.13) 

LaUriX) < ^Ivl^'gix) for X e (6.14) 

Notice that the right hand side in (|6.13p goes to as r ^ since 

r,2 



dr 



<cUynrx\^0 

dr 



Therefore uq is a solution of the homogeneous problem 

uoix, 0) > for X e B'l 

LaUaiX) = for X e Bi\ {{x, 0) : uaix, 0) = 0} 

LaUo{X) < for X e Bi 

For this problem, Almgren's frequency formula applies in the usual way without an error 
correction So, from (j6.1ip . we conclude that uq is homogeneous in -B1/2, and its degree of 
homogeneity is exactly ($u(0) ~ {n + a))/2. Since it is homogeneous, then it can be extended to 
]R"+^ as a global solution of the homogeneous problem. 

The aditional fact that we can choose so that Ur^, \/xUr^ and \y\°'dyUr^ converge uniformly 
comes from the apriori estimates of Proposition 14.31 □ 



Proof of Lemma \6.1[ Two things may happen: 



liminf I ^ ^^^^ ''^^^ 

r^o 1 < +00 second case 
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For the first case, we use Lemma 16.21 to find the blowup profile uq. Then we apply Theorem 
l5.7l to obtain that the degree of homogeneity of mq is 1 + s or at least 2, and thus $„(0) = <&uo(0) = 
n + a + 2{l + s) or $„(0) = $«o(0) > n + a + 4. 

Now we turn to the second case. 

If Fu{rj) < r"+°+"' for a sequence rj 0, then ^ui^j) equals to n + a + 4 for those values of 
rj and <i>„(0) = n + a + 4. 

On the other hand, let us assume that Fu{r) > r"+"+'' for r small. Since we are considering 
the second case, we have 

for some constant C and a sequence rj \ 0. Taking logs in the above inequality 

{n + a + 4) logrj < logF„(rj) < C + {n + a + 4) logrj 

We will show that <f>„(0) = linij^+oo ^ui^j) > n + a + 4 in this case. Suppose otherwise that for 
small rj, $„(rj) < n + a + A — Sq, then we take < r„ << 1 and obtain 

(n + a + 4)(logr„ - logr„i) - C < logF„(r„) - log i^„(r,„) 

< J log Fuir) dr 

< (n + a + 4 - £o)(logr„ - logr„) 

which gives a contradiction when we take logr„ — logr,„ — > +cxd. This finishes the proof. □ 

Proposition 6.3. Let u be as in Lemma \6.1i Assume $(0) = n + a + 2(1 + s). There is a family 
of rotations such that Ur o Ar converges to the unique global profile uq of degree 1 + s (as 
mentioned in Provosition \5.5\) . The convergence is in the sense that 

mW''\B,/2,\yn (6.15) 

uniformly in -B1/2 (6.16) 

uniformly in -B1/2 (6-17) 

uniformly in -B1/2 (6.18) 

Remark 6.4. By a rotation A, naturally what we mean is an orthonormal matrix, i.e. A'^A = Id. 
The set of orthonormal matrices is compact, therefore Proposition 16.31 immediately implies the 
existence of convergent subsequences. 

Proof. In the proof of Lemma 16. li we saw that in the case when liminfr^o < +00 we have 
$(0) > n + a + 4. If $(0) = n + a + 2(1 + s) we must have lim inf = +00. 

Consider = \Hm{B,,,M'^) + Mco{b,,,) + W"^ M\co{b,,,_) + II l2/l''^a^^llc«(Bi/2)- 

From Lemma 16.21 we can find a sequence r^ such that converges to a global profile 

of degree 1 + s. Notice that /^^ |uoP Ij/l^dX — 1, and from Proposition 15.51 uq is unique up to 

rotation. In particular 

lim inf ||mo — u^^ ° ^|| =0 

k — ^00 all rotations A 

where we consider \\v\\ = 11^^11^1(5,/,, |.y|a) + lkllco(i3i/2) + ll^^"llc«(Si/2) + M^'^v'^Wchb,^,)- 

If there was any sequence rk for which ||uo — o A\\ stays away from zero, then by Lemma 
16.21 we can find a subsequence rj-j — > such that converges to a global profile of degree 



Ur O Ar — > Uo 
Ur O Ar — > Uo 
Vx{Ur O Ar) — > V:eMo 

lyl^dyiur o Ar) l^ra^uo 
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1 + s. Notice that J^^ |i*oP |y|°d(T — 1, and from Proposition 15.51 uq is unique up to rotation. In 
particular 



lim inf 

j — >QO all rotations A 



Uq — Ur^. . O A 







arriving to a contradiction. Therefore ||uo — Urj. ° A\\ cannot stay away from zero for any sequence 
Tfc, and the resuh foUows. □ 



Lemma 6.5. Let F be as in (|3.1|) . // 

F{r) < Cj."+'^+2(i+") (6.19) 
for any r < 1, then u{0) = |Vu(0)| 0, and u is C^^" at the origin in the sense that 

\u{X)\ < 

for \X\ < 1/2 and a constant C depending only on C, a and dimension. 
Proof. Consider u+ ~ max(u, 0) and — max(— u, 0). We have 

Lau+ > |yr.g(x) > -C\ynx\ 
Lau- > -\y\''g{x) > -C\ynx\ 

For some r > 0, let uq be the La-harmonic replacement of u+ in Br. We note that 



= LaUo<La[u++C 

\ 2(n+l + a) 

then by comparison principle 

uo>u+- C'r"^ 

By ([QQ]) . we have 

J \uo{X)\' \yrda{X) = J \u+{X)\^ ly^daiX) < Cr^+a+2{i+o.) 

Since jt/l" is an A2 weight, we use the estimates in [H] to conclude 



sup uo < Cr^^" 

Br/2 

Then u^{X) < uo{X) + < Cr^+" for any X E -Br/2- We can use the same reasoning for u~ 
and conclude the proof. □ 

Lemma 6.6. Let be the frequency formula as in (|3.2p . Assume that $tj(r) \ fi as r \ 0. 
Then 

Fu{r) < Cr^ 



for any r < 1 and a constant C depending only on Fu{l) and the constant Co in p.2|) . 

Proof. For simplicity, we omit the subscripts and write <i>(r) and F{r). 

Let F{r) = max(F(r), r"'+''+^). By definition F > F. From the fact that $ is monotone 
nondecreasing 

= $(0) < $(r) = (r + Cor^)-^ logF(r) 

dr 
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then 

— log^^(r)> 



dr r + C'or"^ 

and integrating we obtain 

1 



logF(l)-logF(r) > / J" ^ ds 
J s + Cqs^ 



1 

1 1 



r 

> ~fi\ogr -\og{l - Co)fJ. 

and then log i^(r) < /ilogr + Co/i + logF(l). Exponentiating the two sides of the inequahty we 
obtain 

F{r) < F{r) < Crf" 

where C = F(l)(l + Co)^. □ 
Theorem 6.7. Let u be a function for which ((^ . ([^ . (1^ . (^3)) /loZd and u(0) = 0. Then 

u{X) < C\X\^+'sup\u\ 

Bi 

for some constant C depending only on n, a, and the Lipschitz norm of g. 

Proof. From Lemma 16.11 fi = *i'(0) > n + a + 2(1 + s), so the Theorem is proved by applying 
directly Lemma 16.61 and Lemma 16.51 □ 

Corollary 6.8. If ip G C^'^, the solution u to the problem (jl.ip is in the class C^'^. 

Proof. Using the equivalence between the problem p.ip and p.6m.9p , Theorem 16.71 shows that 
u — (p has a C^'* decay at free boundary points. As shown in [16^, this is enough to prove that 
ueCi'". □ 

Corollary 6.9. If the function u solves the problem i2. 3^2.5]) . then u{x,yo) £ C^''*(i3^^2) ^^2/ 
value < j/o < 1/2. 

Corollary 6.10. // the function u solves the problem il.KMl.f^) . then u(x, yo) G C^''^{B'^^2) for 
any value < yo < 1/2. 

Remark 6.11. As in Remark l4.5[ it is not true that u is C^'* in the whole variable X for general 
values of a. Interestingly enough however, the right decay in y takes place at free boundary points. 
In the variables x and z (as in Remark 14. 5p we would not have the right decay at free boundary 
points, and therefore the optimal regularity would not be uniform in both variables. 



Remark 6.12. The fact that the C^'" estimates for a small happen naturally in the x, z variables 
and the right sharp decay at the free boundary points happens in the x, y variables may somehow 
reflect the respective nondivergent and divergent nature of each of these two results. 
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7 Free boundary regularity 



In this section wc study the regularity of the free boundary. First we show it is a Lipschitz surface 
around nonsingular points. Then we wiU apply the boundary Harnack principle to tangential 
derivatives in order to show that it is a C^'°' surface. This is the same approach that was used in 
[2] and also g]. 

7.1 Lipschitz continuity of the free boundary 

Theorem 7.1. Assume fi — $(0) < n + a + 4. Then there exists a neighborhood of the origin Bp 
and a tangential cone r'(0, e„) C M" x {0} such that, for every r e r'(0, e„) we have DtU > 0. 
In particular, the free boundary F{u) is the graph of a Lipschitz function Xn = fi^i, ■ ■ ■ , a;„_i). 

The theorem will follow by applying the following lemma to a tangential derivative h = D-rUr, 
where Ur is the blow up family that defines the limiting global profiles, for r small. 

Lemma 7.2. Let K be a subset ofR" x {0}. Assume h is a continuous function with the following 



4- h> —uj{a) for \y\ < a, where uo is the modulus of continuity of h. 

There exists cto = cro(?i, a, cq, lo) and 70 = "io{ri, a, cq, uj) such that, if a < aa, 7 < 70 then h > 
in Bi/2- 

Proof. Suppose Xq = {xo,yo) € -B1/2 and h{Xo) = 0. Let 



properties: 



1- \Lah\<-i\y\ 



a 



2. h>0 for \y\><J>0,h = on A. 



3. h>cn>0 for \y\ > 





and 



P{x,y) = \x- XoP 



n 




a+1 



Observe that LaP — 0. Define 



v{X) = h{X) + SP{X) 



7 




2(a+ 1) 



Then: 



v{Xo) = h{Xo) + SP{X„) - - 




LaV = Lah + SLaP ~ i\yT < < 
v{X) > 



outside A 
on A 



Thus, V must have a minimum on dQ. On 
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we have 



S 7 

v>cq — > 

16n 32n 



if S and 7 are small depending on cq and n. 
On 



, 1 , , VT 

F - a^ol = o ^ \y\ ^ 



one has 



if 7 is small compared to 6. 
Finally, on 



3 - i«i - 8„ 



„>^( 1- — ) ^>0 

- 9 64nJ 128n2 



a^ol = ^ |y| < cr 



one has 



if a is small, depending on S and n. 

Hence, u > on dQ and we have a contradiction. Therefore /i > in J5i/2- D 

Theorem 17.11 follows applying the Lemma [7.21 to a tangential derivative h = DrUr, where Ur 
is the blow up family that defines the limiting global profiles, for r small (also to keep small the 
constant 7 appearing in the lemma). This is a standard procedure by now that mimics the proof 
of Lipschitz regularity of the free boundary in the classical obstacle problem [4] . Below we include 
the proof for completeness. 

Notice that by the properties of the limiting profile we can actually choose the opening 9 of 
the cone as large as we like. 

Proof of Theorem \7J\ We have /i = $(0) < n + a + 4. By Lemma 16. 1[ /i = n + a + 2(1 + s). 
Moreover, the blowup sequence Ur converges (up to subsequence) to the global profile uq of degree 
1 + s for which the free boundary is fiat. 

Let us assume that e„ is the normal to the free boundary of uq- We know by Proposition 15.51 
that ^ 

dnUo = C (^y^xl + y2 - Xn^ 

For some 6*0 > 0, let v be any direction orthogonal to y and x„ such that \v\ < Oq. From 
Theorem 15. 71 we have that uq is constant in the direction v, therefore if r = e,i + v, OtUq — dnUQ. 

We know by Proposition 16.31 that Vj^u^ converges uniformly to V^uo (up to subsequence). 
Therefore, for any 60, there is an r for which \dvUr — drUo\ < So for any r constructed as above 



(^0 depends on 9q). We also have that drUo — 9„Wo ~ c y\/ x"^ + y'^ — x„j . If we differentiate 
(|6. 12116. 14|l , we obtain the equations for Ur 

La{drUr{X)) = '^\y\-rdrg{rx) < Crlyj" for X&Bi\ {(x,0) : Ur{x, 0) = 0}, (7.1) 

and its right hand side goes to zero as r — > 0. Therefore for a small enough r, dfUr will satisfy 
the hypothesis of Lemma 17.21 Thus drUj. will be nonnegative in Bi/2- This implies that for any 
point X in the free boundary of u such that |x| < r/2, then x + Xt stays on one side of the free 
boundary for any positive value of A, and thus the free boundary is Lipschitz. □ 
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7.2 Nondegeneracy properties 

Assuming as before that the origin is a free boundary point, we cannot assure any minimal growth 
of u around without additional assumptions. It is easy to write down examples for the case 
a = in two dimensions. Even for homogeneous global solutions vanishing on half of M" we have 
the following family of solutions of the thin obstacle problem expressed in polar cordinates. 

u = r^''+'^ sm{{2k + 1)9) 

which grow arbitrarily slowly away from the origin. 

The nondegenerate point that we will consider are those such that $(0) = {n + a) + 2(1 + s). 
At those points we know the exact asymptotic picture from Proposition [^31 For some small r, Ur 
will be close to the blowup profile, and we will be able to apply the following lemmas. 

Lemma 7.3. There are positive numbers Eq — eo{a,n) and Sq — (j^)^^^'^ so that the following 
is true: 

Let V be a function such that 

Lav{X) < eo for X e x (0, 6o) (7.2) 

v{X) > forX£ Bl X (0, 5o) (7.3) 

v{x,5^)>^ forx&Bl (7.4) 

Then v{x,y) > Cly]^" in B*/2 ^ i^^So]. 

Proof. This lemma refines the result of Lemma l7.2l and therefore we can expect the proof to follow 
more or less the same pattern. 

Let Xq = (xq. Ho) G ^ ['-'' "^o]- We compare v with the function 

-(x,y) = (l + f)y^-^^^+j/- 
2 n 

We observe that inside the set x (0, 6o), LaW = Sq > LaV. We must check at the boundary 
now. For y = 0, w{x,0) < < v{x,y). For y ~ 5f), 

wix, So) < (1 + §)<52 + < 3(<5o)2^ < < vix, So). 
2 All 

And finally, where \x — xo\ = 1/2, w{x, y) < 3(5q* ^ ^ < < v{x, y). Therefore by comparison we 
have that w < v in the whole set B'l x [0, Sq]- In particular 

v{xo,yo) > w{xo,yo) > 

□ 

Corollary 7.4. Let u be a solution of i2.^2.5]) . Moreover, let us assume that the function g in 
(|2.4p is small in the sense that \g{x)\, \\I xg{x)\ < eo- Let uq be a global profile as in Proposition 
15.51 such that |Va;U — Va;Uo| < eo- Then if eo is small enough, there is a constant c depending only 
on n and a such that we have the following lower bound for the tangential derivatives Ur 

UriX) > cdist(X, A)2* 

for any X G -Bi/g and r being a tangential unit vector such that |t — e„| < 1/2. 
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Proof. From Lemma 17.21 and (|7.ip . we know that Ur will be positive in i?i/2- 
Applying Lemma l7.3l to an appropiate multiple of Ur we immediately get that 



Ur{x,y) > c\y 



2s 



for every {x,y) G i?i/4- Now, let X — {x,y) e Si/g and d = dist(X, A). Consider the ball 
B^I2{X), then for the point T = [xt, yr) at the top of this ball we have yr > rf/2. Therefore 

Ur{T) > Ccf' 

By Harnack inequality fProposition l2.2| ) we know that Ut{X) and Ut(T) are comparable since 
Ur is a solution of the equation in the ball Bd{X). Therefore also 

Ur{X) > cd?" 

□ 

7.3 C^'" regularity of the free boundary 

As in [4] or [2], we intend to prove a C^'" regularity estimate for the free boundary applying the 
boundary Harnack inequality to two positive tangential derivatives Ur- We must be careful in 
order to handle the right hand side in (|2.2ll2.5p . 

The first part in boundary Harnack (the Carleson estimate) can be proved even with a right 
hand side if we have a nondegeneracy condition like the one given by CoroUarv 17.41 

We use the notations from [2] 

• Z) C Si is a La - NT A domain. 

• = dD C] Bi. If Q e il, Ar{Q) is a point such that S,,,.(A((5)) C Br{Q) n D. In particular 

\Ar{Q)-Q\^r. 

Notice that the definition of NT A domain given in [5] differs slightly from the one used in [TT] . 
A uniform capacity condition is assumed for the complement, instead of containing a full ball. 

Lemma 7.5. Let D d Bi be a La — NT A domain. Suppose u is a positive function in D, 
continuously vanishing on f2. Assume that: 

\Lau\ < C\yr (7.5) 
u{X) > cd^ for some f3, < (3 < 2, where dx = dist(A", fl) (nondegeneracy) (7.6) 

The for every Q ^ n B1/2 and r small: 

snp u<C{n, a, D)u{Ar{Q)) (7.7) 
BAQ) 

Proof. Let u* be the La-harmonic replacement of u in B2r{Q) H D, r small. Then by thorem 6 in 
[2] adapted to the operator La (the proof is identical using the results in [12]) we get 

u*{X) < Cu*{AriQ)) in S,.(g) n D (7.8) 

On the other hand, 

u* (X) + C{\X ^ Q\^ ^ r^) < uiX) on d{D n B2r) 

La {u* (X) + C(| A - Q|2 - r^)) = C\y\^ > Lau{X) in n 
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then by comparison principle u* ~ u < in D D i?2r- Similarly we obtain the other inequality 
and \uiX)-u*{X)\ < Cr^. 
From 

u{X) < CiuiAr{Q)) + Cr^). 
From the nondegeneracy property, u{Ar{Q)) > cr^^, and since /? < 2, (|7.7p follows. □ 

We will not prove the second part of Boundary Harnack in so much generality. Instead, we 
notice that for our particular case, we can obtain it from Lemma 17.21 

Lemma 7.6. Let D = Bi — A, where A is a subset o/R" x {0} such that there is a constant k, so 
that the La-capacity of Br(X) Ci A is at least nr"^"-^^ for every X A and r < 1/2 ( this implies 
that D is an NTA-domain as in 

Let u,v be two positive functions in D, as in Lemma \7.5\ Let us also assume that u and v are 
symmetric in y. Then there is a constant c (depending on n, a, and the constant k above) such 
that 

u{X) u{en+i) , . 

-TFT < c-^ r in Bi/2 (7.9) 

v{X) v{en+i) 

Moreover, the ratio u/v is a C" function in -B1/2 (uniformly up to the boundary A). 

Proof. We divide both functions by their values at e„+i, so that we can assume u(e„+i) = 
v{en+i) = 1. 

Using Lemma 17.51 and Harnack inequality, for any 6 > we have 

u{X) < C for aU X e B3/4 

v{X) > c for aU X e B3/4 n {y > 6} 

which means that for a constant s > small enough, v — su satisfies the assumptions of Lemma 
17.21 Therefore i; — su is positive in B1/2, or in other words u{X)/v{X) < s. 

By a standard iterative argument (|7.9p implies that u/v is C" up to the boundary for an 
appropriately small a. □ 

The C^'" regularity of the free boundary follows by applying Lemma 17.61 to two positive 
tangential derivatives. We need the obstacle to be C^-^. 

Theorem 7.7. Let u be a solution of Il2.0i2.5\) . Assume is in its free boundary and $(0) < 
n + a + 4. Then the free boundary is a C^'°' {n — 1)- dimensional surface around 0. 

Proof. The proof follows by a standard argument in the subject [4], [^. 

From Theorem 17.11 the free boundary is Lipschitz. Moreover there is a cone of tangential 
directions r, so that Ut is positive in a neighborhood of the origin. Let us assume that e„ is an 
axis of that cone. We express the tangential derivatives of the level sets of m as a quotient u-r/ue„ 
and use Lemma l7.6l to show that all level sets up to the free boundary are uniformly C^'°' surfaces 
in the tangential direction. 

We must show that we have all the necessary conditions to apply Lemma [7.61 First of all, we 
consider a rescaled function Ur so that Ur is close enough to the global profile uq and the right 
hand side in (|6.13p is small. 

If T is in the cone of directions from the free boundary where the directional derivatives are 
positive, then Ur is positive. Otherwise we express 

so we consider e„ + t/4 instead, and reduce the problem for those r where Ut is positive. 
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From Corollary 17.41 we have the nondegeneracy condition for both Ue„ and . In order to 
apply Lemma [7.61 we are only left to show that the contact set A = {(x, 0) : u{x,Q) — 0} has a 
uniform La-capacity. But that is automatic since the free boundary is Lipschitz. □ 

Appendix: Frequency formula. 

In this appendix we prove Theorem 13. II We recall the notation from (|3.ip 



where u is a function that solves (|2?2l) . (fO]) . ([2^ . (|2?5l) , and u(0) = 0. 

As it was mentioned before, in the case of zero right hand side in ()2.4|) . the frequency formula 
takes the simplest form and was proved in [7j. In general we will show that a constant Co exists 
so that 

$(r) = (r + Cor2)4-logmax(F„(r),r"+'^+'^) 
dr 

is monotone increasing. 

In order to prove that $ is monotone, we will need the following lemma 

Lemma 7.8. The following identity holds for any r < 1. 

r J {\ur\^ - \u,f) lyrdd = I ((n + a - 1) \Wuf - 2{X, Wu)g{x)) \y\^dX (7.10) 

where Ur is the gradient in the tangential direction to Sr and u^, is the derivative in the normal 
direction. 

Proof. The Lemma is obtained by applying the divergence theorem to the vector field 

-X- \y\''{X,Vu)Vu 



We have to notice that 

,„,(„.^._,»r<.-,v..,v„)...ii±|^,v„p-<x,v„,..„ 

Since (A", Vu) is a continuous function on B* that vanishes on A = {u = 0}, we have that 
{X, \Ju)LaU has no sigular part and coincides with {X, Vu)g{x). □ 



We now prove Theorem 13. II 
Proof of Theorem \3.1\ Notice that log i^(r) is differentiable for r > since 

— logi^(r) - ^'^''^ - ^/s,-^^'- 1^1°^'^ ^ n + a 
dr F(r) \y\°-d(7 r 

On the other hand logr"+"+'* clearly has "^"''"^ as derivative. When we take the maximum 
in (|3.2p it may happen that we get a nondifferentiable function, but in any case the jump in 
the derivative will be in the positive direction. In order to prove that $ is monotone, we can 
concentrate in each of the two values for the maximum separately. 
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In the case $(r) = (r + Cor^)-^ logr"+°+'*, we have $'(r) = (1 + 2Cor)(n + a + 4). Therefore 
$ is monotone in this case. 

When $(r) = {r + Cor'^)-^ log F{r) we need to work a httle more. We will concentrate in this 
case now. Notice that this case happens only if F{r) > 7-"+'^+"'. 

We expand the expression for ^{r) 

The second term {n+a){l + Cor) is clearly increasing, so it would suffice if we can show that the 
first term is nondecreasing. We compute its logarithmic derivate and show that it is nonnegative, 



■log (r + Cor^)- 



dr \ Jg \y\°'da j 

log(r) + log(l + Cor) + log j uu, \y\''da - log j v? jy^da (7.11) 



_d 

dr 



Sr Sr 

Co -^/s,""i^ \g2uuy\y\''da n + a 



r 1 + Cor uui, |y|°dcr lyl^dcr r 

To estimate uuy \y\°'da we notice that 

y uu^ |y|"dCT = y divdyl'^uVu) dX 



|y|"|Vu|^ + uLaudX 



Thus 



d 
dr 



(7.12) 



y uu,. |yrdo- = J \y\°-\Vu\' + uLaU da 

> I |Vu|2 |y|"da-Cr^ (j \u\' lylMaj (7.13) 

> J |Vu|2 |y|M(T-Cr'^F(r)^/2 

Sr 

Now we use Lemma FTSl to estimate iVwp |y|°d(T. 
|Vw|'|yrda = 2 /" |u,ny|"d(7 + i / ((n + a - 1) |Vu|' - 2(X, Vu)5(a;)) jy^dX 



Br 

= 2 J \u^f |y|°dcr + ^^'^^^ J uu^ |y|°dcr - ^ ^ ii^ + a~-l)u + 2{X, Vu))g{x) ly^dX 

Sr Sr Br 



28 



Recall that \g{x)\ < C\x\. Putting the above estimate back in (I7.13P we obtain 

Sf Sf S r 



where 



F{r) = J u{Xf lyrda 
G(r) = j u{Xf l^rdX 

B,. 

H{r) = j |Vw(X)|2 \y\-dX 



Putting it all together in (|7.1ip . we write log<i>(r) = R{r) + S{r), with 

^ 2/gJz..p \y\-da _ /^^ 2uu, |y|"da ^ ^ 
~ uu^ lyl-^dCT /^^ 1*2 \y\''da " 

and 

1 + Cor UM^ lyhdcr 

In order to estimate S properly, we must find bounds for _F, G and H . Recall that we are in 
the case when F{r) > From Lemma [2. 131 we know 

^n+a+i < p(^^^ < CrH{r) + Cr^+'^+" 
and also, integrating the inequality in r we have 

G{r) < Cr^H{r) + Cr''+^+" 
This means that, for small enough tq, if r < tq, 

F{r) < CrH{r) 
G{r) < Gr^H{r) 

and then 



1 + Cor H{r)~-r''^y/H{r) 
But since F(r) > also > We thus have, for r < ro, 

^ ^ - 1 + Cor VW) 

-1 + Cor 

>^^~C 
-1 + Cor 
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Which is positive if we choose Co large enough. This shows that $'(r) > for smaU r. If 
$ is not monotone for larger values of r in (0,1), we just take Co larger enough so that $ is 
nondecreasing for < r < 1. □ 
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